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ABSTRACT 

We present high speed photometric observations of the eclipsing dwarf nova IP Peg 
taken with the triple-beam camera ULTRACAM mounted on the William Herschel 
Telescope. The primary eclipse in this system was observed twice in 2004, and then a 
further sixteen times over a three week period in 2005. Our observations were simul- 
taneous in the Sloan u' , g 1 and r' bands. By phase-folding and averaging our data we 
make the first significant detection of the white dwarf ingress in this system and find 
the phase width <f> of the white dwarf eclipse to be 0.0935 ±0.0003, significantly higher 
than the previous best value of 0.0863 < (j> < 0.0918. The mass ratio is found to be 
q = M2/M1 = 0.48 ± 0.01, consistent with previous measurements, but we find the 
inclination to be 83.8 ± 0.5 deg, significantly higher than previously reported. We find 
the radius of the white dwarf to be 0.0063 ± O.OOO3i?0, implying a white dwarf mass 
of 1.16 ± 0.02M Q . The donor mass is 0.55 ± 0.02M Q . The white dwarf temperature 
is more difficult to determine, since the white dwarf is seen to vary significantly in 
flux, even between consecutive eclipses. This is seen particularly in the u'-band, and is 
probably the result of absorption by disc material. Our best estimate of the tempera- 
ture is 10, 000 - 15, 000K, which is much lower than would be expected for a CV with 
this period, and implies a mean accretion rate of < 5 x 1O _11 M0 yr _1 , more than 40 
times lower than the expected rate. 
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1 INTRODUCTION 



Cataclysmic variable stars (CVs: Wa rnerlll995l ) provide ex- 
amples of white dwarfs accreting from low mass compan- 
ions at rates of ~10~ n - 10" 9 M Q /yr. The mass transfer 
in CVs and similar semi-detached binary systems is driven 
by angular momentum loss. For long orbital period (~3 h 
or greater) systems the angular momentum loss is thought 
to be driven by magnetically-coupled stellar winds. For the 
shorter period (< 3 h) systems losses due to gravitational 
radiation are thought to dominate the mass transfer pro- 
cess. This is observationally supported to some extent: the 
shorter period population is dominated by lower M dwarf 
novae, whereas the peak of the populati on of high M 'steady 
state' CVs is in the 3 to 4h period range (IShafterll992T ). How- 
ever, there is no one-to-one relationship. The distribution of 
steady-state, bright CVs tails off sharply at longer periods 
jRitter fe Kolbll2003l h and a number of dwarf novae are ob- 
served to have very long periods. It is generally thought that 
all non-magnetic CVs are essentially the same, with only a 
change in M required to move between the different sub- 



classes. While it may be the case that the high- and low-M 
systems at a given orbital period are different populations 
with a different evolutionary history, a more likely expla- 
nation is that the instantaneous M we measure is a poor 
indicator of the mean M, so that even high mean M sys- 
tems spend some fraction of their time as slowly accreting 
dwarf novae. 

One independent means of inferring a longer-term (10 3 
- 10 yr) value of M is through measurement of the white 
dwarf temperature. The temperature is a good tracer of 
the long-term M since it is determined not by the ac- 
cretion heating, but by the compress ion of the underlying 
white dwarf by the accreted mat ter (|Townslev fc Bildstenl 
2003;; ITownslev fc Gansickd | 20091 ) . The existing measure- 
ments agree well with the high -M/long period, low-M /shor t 
period consensus (figure 5 of ITownslev fc Gansickej |2009). 
However, these data are likely to be biased. These white 
dwarf temperatures are spectroscopicalfy determined, and 
the white dwarf must dominate at UV wavelengths in order 
for the temperature to be determined. Long period CVs are 
generally brighter than shorter period systems due to their 
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higher accretion rates and more luminous donor stars, and 
so the white dwarfs in long period systems must be hotter 
in order to be measurable. An alternative strategy which 
reduces this bias is to make fast, multi wavelength photo- 
metric observations of CVs where the inclination is high 
enough for the white dwarf to be eclipsed. With sufficient 
time resolution the ingress and egress of the white dwarf 
can be identified independently of the other components, 
and hence the white dwarf colours can be determined. 

An ideal system for a study of this nature is IP Pegasi 
(IP Pe g hereafter). Thi s dwar f nova was found to be eclips- 
ing by iGoranskii et al] (|l985l ) and has a period of 3.8h, in 
the 'lo ng period' tail of the dwar f nova distribution. Fol- 
lowing |^^nsleiTBiidste3 (|2003l ) this would imply a high 
mean M and a correspondingly high white dwarf tempera- 
ture. However, the temperature of the white dwarf has yet 
to be determined with precision. Previous studies of IP Peg 
have shown that the white dwarf ingress feature is obscured 
by the ingress of the bright spot on the accretion disc, and it 
has been claimed that the white d warf egress feature varies 
considerably in size and length (|Wood fc Crawford! 1 19861 
I Wood et all 1 19891 ). The difficulty in detecting the white 
dwarf throws the system parameters of IP Peg into doubt. In 
an attempt to resolve this and to measure the temperature of 
the white dwarf we took a total of eighteen separate observa- 
tions of the IP Peg eclipse with the high speed CCD camera 
ULTRACAM mounted on the William Herschel Telescope 
(WHT). In this paper we present these data, and our sub- 
sequent parameter determinations. 



2 OBSERVATIONS 

The high speed CCD camera ULTRACAM l|Dhillon et all 
120071 ) was mounted on the 4.2m William Herschel Tele- 
scope (WHT) in August 2005. Observations were made with 
this instrument between 9-15 August, and 25 August - 1 
September. Over these periods, 16 separate observations of 
the IP Peg eclipse were made. In addition, two observations 
of the eclipse were made in August 2004. ULTRACAM is 
a triple beam camera and all observations were made us- 
ing the SDSS u , g' and r' filters, except for the first 2004 
observation, in which the i' filter was used in place of r', 
for scheduling reasons. Most observations were taken with 
an exposure time of 1.5s, which gave a sufficient count rate 
in all three bands and adequately sampled the white dwarf 
egress, which lasts ~30s. The dead time between exposures 
for ULTRACAM is ~25ms. The data were unbinned and 
the CCD was windowed in order to achieve this exposure 
time. This exposure time was varied by up to Is in order 
to account for changing conditions. A complete log of the 
observations is given in Table [1] 

All of these data were reduced with aperture photom- 
etry using the ULTRACAM pipeline software, with de- 
biassing, flatfielding and sky background subtraction per- 
formed in the standard way. The source flux was deter- 
mined using a variable aperture (whereby the radius of the 
aperture is scaled according to the FWHM by a factor of 
1.7). Variation in observing conditions were accounted for 
by dividing the source lightcurve by the lightcurve of a 
nearby comparison star. We use the star at a = 23/i23m06s, 
5 = +18°24'40" as the comparison. We checked the stability 



of this star against other stars in the field, and calculated 
apparent magnitudes of 14.45, 13.20 and 12.62 in u' , g' and 
r' respectively. 

Our short eclipse observations were interleaved with 
longer observations for other programmes. We used the field 
stars in these long datasets to determine atmospheric ab- 
sorption coefficients in the u' , g' and r' bands, and subse- 
quently determined the absolute fl ux of our targets us ing ob- 
servations of standard stars (from S mith et al.ll2002l ) taken 
in evening twilight. We use this calibration for our deter- 
minations of the apparent magnitudes of the two sources, 
although we present all light curves in flux units determined 
using the conversion given in ISmith et all (|2002l ). Using our 
absorption coefficients we extrapolate all fluxes to an air- 
mass of 0. For all data we convert the MJD times to the 
bary centric dynamical timescale, applying light travel times 
to the solar system barycentre. 



3 LIGHTCURVES 

In Figure [1] the eighteen eclipses that were observed in 2004 
and 2005 are plotted. In some cases we only have a partial 
lightcurve, due to either poor weather or scheduling con- 
straints. We comment on weather conditions in Table [1] In 
particular, note that the observations made on 25 August 
and later (eclipses 9 to 18) suffered from high levels of Sa- 
haran dust in the atmosphere. This dust has a grey absorp- 
tion curve but the loss of flux is particularly apparent in the 
w'-band lightcurves. 

These lightcurves show a single ingress feature, followed 
by separate egresses of the white dwarf, and later, the bright 
spot. While the egresses of these two components are dis- 
tinct, it is impossible in these plots to disentangle the two in- 
gresses. There is considerable accretion-driven flickering out 
of eclipse, most of which is due to the bright spot, although 
in some cases there is a degree of variation immediately af- 
ter the white dwarf egress. There is noticeable variation in 
the strength of the white dwarf egress feature. This will be 
addressed in detail in Section ^. 21 The 2004 observations are 
at a different point in the outburst cycle to the 2005 obser- 
vations, and so the shape of eclipses 1 and 2 are different to 
the rest, due to the accretion disc flux. In addition, in eclipse 
1 there is a 'hump' between the white dwarf and bright spot 
egresses, which may be due to flickering. Finally, it can be 
seen in eclipse number 3 that there is a brief brightening 
in the lightcurve just prior to the white dwarf egress. This 
brightening is most likely a flare event from the donor star, 
and makes it impossible to model the egress for that eclipse 
accurately. 



4 ECLIPSE MODELLING 

In this section we report system parameters determined by 
fitting the data with the model described in Appendix [X] 
This was a three stage process. We first obtained an ini- 
tial fit to each individual eclip se using th e simplex and 
Levenberg-Marquardt methods (Press 2002). These initial 
fits were used to obtain eclipse timings, and the resultant 
ephemeris was used to phase-fold all of the complete eclipse 
lightcurves. This combined lightcurve was used to determine 
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Table 1. Log of the observations. Each observation is a single eclipse of IP Peg. Unless otherwise stated, the Sloan u' , g' and r' filters 
were used. We give the phase range centred on the white dwarf egress. 

Eclipse Date UT Exposure Phase 



number 


Cycle 


start 


start 


end 


time (s) 


ran 


ge 






Comments 


1 


-2200 


29 Aug 2004 


01 


10 


06 


05 


1.5 


-0 


60 





70 


Clear, seeing 0.6 - 0.8" . 
i'-band filter used in place of r' . 


2 


-2194 


30 Aug 


02 


04 


02 


24 


1.5 


-0 


06 


- 


02 


Clear, seeing 0.6 - 0.8". 

Partial lightcurve: White dwarf egress only. 


3 


-1 


12 Aug 2005 


00 


41 


01 


17 


1.5 


-0 


09 


- 


07 


Clear, seeing 0.6 - 1.0". 

M-dwarf flare just before white dwarf egress. 


4 





12 Aug 


04 


17 


05 


47 


2.5 


-0 


11 


- 


21 


Mostly clear, seeing 0.6". 


5 


6 


13 Aug 


03 


06 


04 


00 


1.5 


-0 


10 


- 


10 


Clear, seeing 0.6". More readout noise than 
usual in the r'-band due to a CCD problem. 


6 


18 


15 Aug 


00 


33 


01 


34 


1.5 


-0 


16 


- 


10 


Clear, seeing 0.6". 


7 


19 


15 Aug 


04 


20 


05 


31 


1.5 


-0 


03 


- 


11 


Clear, seeing 0.6". 

Partial lightcurve: white dwarf egress only. 


8 


21 


15 Aug 


23 


17 


00 


17 


1.5 


-0 


17 


- 


08 


Clear, seeing 0.6". 


9 


87 


25 Aug 


22 


28 


23 


27 


1.5 


-0 


08 


- 


08 


Partial lightcurve due to cloud. 
Seeing 0.6 - 1.0". Dust. 


10 


88 


26 Aug 


01 


47 


03 


18 


1.5 


-0 


18 


- 


09 


Partial lightcurve due to cloud. 
Seeing 0.6 - 1.0". Dust. 


11 


101 


28 Aug 


03 


30 


01 


35 


0.9 


-0 


21 


- 


07 


Seeing 0.5". Dust. 


12 


106 


28 Aug 


22 


43 


23 


38 


1.2 


-0 


15 


- 


08 


Seeing 1 —2". Some dust. 


13 


107 


29 Aug 


02 


24 


03 


21 


1.2 


-0 


17 


- 


07 


Seeing 1 —2". Some dust. 


14 


112 


29 Aug 


21 


47 


22 


18 


2.5 


-0 


07 


- 


06 


Seeing 1 —3". Dust. 

Partial lightcurve - no ingress. 


15 


113 


30 Aug 


01 


18 


02 


10 


1.5 


-0 


11 


- 


08 


Seeing 1 — 2". Dust. 


16 


119 


31 Aug 


00 


02 


00 


57 


0.9 


-0 


16 


- 


08 


Seeing ~0.8". Dust. 


17 


120 


31 Aug 


03 


50 


04 


43 


0.9 


-0 


16 


- 


07 


Seeing 0.5 - 0.6". Much dust. 


18 


133 


02 Sept 


04 


25 


06 


11 


1.2 - 2.0 


-0 


36 


- 


08 


Seeing 1 — 2". Dust. 



the mass ratio q and inclination i of the system. We use a 
Markov Chain Monte Carlo (MCMC) algorithm for minimi- 
sation and determination of uncertainties. These parame- 
ters were used as a basis for MCMC fits of the individual 
lightcurves, so as to determine parameters which would be 
expected to vary over the course of our observations, such 
as the accretion disc radius. Since the white dwarf egress 
feature appears to vary in height, and since there have been 
reports of variable egress duration l|Wood fc Crawfordlfl986l ; 
I Wood et al.ll 19891) , the white dwarf flux and radius were de- 
termined individually for each eclipse, as well as more refined 
eclipse timings which we used to determine a new ephemeris 
for this system. 

4.1 Fitting the phase-folded data 

We selected the eleven lightcurves in 2005 where the entire 
eclipse is covered (numbers 4, 5, 6, 8, 11, 12, 13, 15, 16, 17 
and 18), and phase- folded them using the ephemeris given 
in Section 14.31 The combined and binned lightcurves are 
plotted in Figure [2j using all eleven complete eclipses. Each 
bin is ~2.7s in length. With the flickering reduced, a white 
dwarf ingress feature is clearly apparent, distinct from and 
preceding the ingress of the bright spot. We measured the 
phase width of the white dwarf eclipse to be 0.0935 ±0.0003, 
which is great er than the value of 0.0863 suggested as being 
most likely by I Wood fc Crawford! (|l986h and is also greater 
than the upper limit of 0.0918 given in the same paper. This 
implies either a higher orbital inclination or a higher mass 
ratio or some combination of the two. 

It might be thought that this feature could be due to 



nickering. However, in all three bands it has a height and 
duration that is consistent with it being the white dwarf 
ingress. Additionally, by changing the number of lightcurves 
used to create the phase-folded plot, we find that the ex- 
istence and size of this feature in the combined dataset is 
not due to any one individual lightcurve. This is not true of 
other significant variations in the lightcurve, such as the two 
peaks immediately preceding the ingress at a phase range of 
—0.11 < <f) < —0.10, which are definitely due to nickering. 
We are therefore confident that we have detected the ingress 
feature. 

A model consisting of a white dwarf primary, a Roche 
lobe filling secondary and accretion disc and bright spot 
components was fitted to these data using the MCMC 
method. A complete description of this model is given in 
Appendix [X] In Table lATI we list the parameter determina- 
tions we make with this method, although some of these are 
refit on an eclipse-by-eclipse basis in Section \4. 2 1 The uncer- 
tainties listed in this table are the formal errors, scaled so as 
to give a reduced x 2 of 1. These uncertainties are underes- 
timates, as shown by the inconsistencies between the three 
photometric bands. In Section f5.1l we will discuss this issue 
further, and calculate more realistic uncertainties which ac- 
count for the systematics. The best model fits are plotted 
in Figure [2] The fits are generally good, although there is 
some variation in the residuals out-of-eclipse due to flick- 
ering which has not been completely averaged out in the 
combined lightcurve. Additionally, our fit to the bright spot 
ingress and egress is not perfect; we discuss this in Appendix 

m 

One additional parameter which our results are sen- 
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Figure 1. Lightcurves of IP Peg eclipses, observed in August 2004 and August 2005. We label the eclipses according to the numbers 
given in the observation log (Table [TJ. All data were collected simultaneously in the u'- g 1 - and r'-bands (blue, green and red lines, 
respectively) except for eclipse 1, in which an i' filter was used in place of r'. Each lightcurve in this plot is offset from the previous by 
1.5mjy. 



sitive to is the limb darkening coefficient for the white 
dwarf. An initial fit to the data was made with the 
Levenberg-Marquardt method in order to obtain approx- 
imate white dwarf parameters, and these determinations 
were used with the white dwarf m odel atmospheres of 
iGansicke. Beuermann fc de Martinol (|l995T l to determine 



the coefficients in the three ULTRACAM bands. We found 
linear coefficients of 0.352, 0.253, 0.216 in u' , g' and r' to 
be appropriate for the white dwarf in IP Peg. These values 
were used in all our MCMC fits. 
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Figure 2. Top panel: phase-folded and binned lightcurves (top, r'\ middle, g'\ bottom, «'), comprised of the eleven complete eclipse 
observations made in 2004 and 2005. We plot the average flux in mjy against the binary phase, where a phase of corresponds to the 
mid-egress of the white dwarf. We plot the datapoints with uncertainties in black, and the best model fits to these data in red, green 
and blue (for r', g 1 and u', respectively). The four vertical lines indicate the start of (from left to right) the white dwarf ingress, the 
bright spot ingress, the white dwarf egress and the bright spot egress. Middle panel: residuals in the three bands (top, r'\ middle, g'; 
bottom, «'). Bottom panel: the three components of the g'-band model plotted separately, showing the relative strengths of the bright 
spot, accretion disc and white dwarf. These three lines have been offset for clarity. 



4.2 Fitting the individual eclipses 

For some parameters the results from the phase-folded 
lightcurves are insufficient. Some variations in the white 
dwarf egress feature are clearly apparent in the individual 
eclipses, so it was necessary to fit these separately in or- 
der to obtain a proper determination of the white dwarf 
radius and fluxes. Obviously this was also necessary in or- 
der to determine the eclipse timings, and some variation 
in the accretion disc radius over the course of our observa- 
tions would also be expected. We used our model fit to the 
phase-folded lightcurve as a basis for fitting the individual 
eclipses separately. It is clear from Figure [T] that the ampli- 
tude of the flickering is much larger than the white dwarf 
eclipse features in these individual lightcurves, so we fixed 
most of the parameters, including the mass ratio and incli- 
nation angle, to the values determined from the phase-folded 
lightcurve. We allowed the white dwarf radius and the time 
of the white dwarf mid-eclipse to vary. We also wanted to 
track variations in the disc radius, but in these individual 
lightcurves the outer disc radius Rdisc is poorly constrained. 



We therefore set Rdisc = Rspot and fitted R 3po t , where R apo t 
is the distance between the white dwarf and the maximum 
of the surface brightness of the bright spot. This parame- 
ter is well constrained by the data, since the bright spot 
is so strong in flux. Note that in the previous section (the 
phase-folded lightcurves) the data are of sufficient quality 
that we could fit Rdisc and R spo t separately, and we found 
Rspot to track the accretion disc radius closely. To allow for 
a degree of flexibility in modelling a variable bright spot, we 
also allow I (the scale length of the bright spot) and f c (the 
fraction of bright spot flux that is taken to be constant) to 
vary. The remaining bright spot parameters were fixed to 
the values determined from the phase-folded lightcurve. Fi- 
nally the surface brightness of the accretion disc is modelled 
as a power law, and the exponent of this power law is also 
allowed to vary. A complete description of all parameters is 
given in Appendix [X] Each individual lightcurve was fitted 
in the same way as the phase-folded lightcurve (as described 
in Section I4T) . 

The results for the p'-band fits (for which the signal-to- 
noise is highest) are listed in Table [2] We report the white 
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dwarf radii and timings for seventeen of the eighteen eclipses, 
omitting the eclipse with the flare near the egress feature. 
There is a considerable degree of scatter in our model de- 
terminations of the white dwarf radius. However, given the 
uncertainties, the variations are not significant. The timings 
of the white dwarf mid-eclipse were converted to timings 
of the mid-point of the egress feature so as to be consis- 
tent with previous authors, using a phase width of 0.0935. 
We list also our determinations of the accretion disc radius 
Rdisc scaled by the binary separation a. This parameter is 
given for fifteen of the eighteen eclipses, omitting eclipses 
which are not complete enough for the disc radius to be de- 
termined with any precision. Also listed in Table [2] are the 
white dwarf fluxes in all three bands. These were calculated 
by measuring the height of the white dwarf egress feature 
in the model fi ts and then conver ted back to a magnitude 
scale following ISmith et all <|2002h . In Figure [3] the u'-, g'- 
and r'-band lightcurves of the region around the white dwarf 
egress are plotted, with the accretion disc component of the 
emission subtracted from the data. This plot shows clearly 
that the height of the egress feature varies significantly from 
eclipse to eclipse (for example, compare egresses 15 and 16, 
taken on consecutive nights). This variability is apparent in 
all three bands, although in r' the egress feature is very small 
and in v! the signal-to-noise ratio is much lower, making it 
harder to determine the size of the egress accurately. 



4.3 Eclipse timings 

We used the timings of the white dwarf mid-egress 
given in Table [5] in order to calculate an updated 
orbital ephemeris for IP Peg. These data were comple- 
mented with eclipse timings g iven b y previous authors 
dGoranskii Lvutvi fc Shugarovl Il985l ; IWood et al.1 1 19891 ; 
IWolf et al.lll993l ) All times were converted to the barycen- 
tric dynamical timescale, corrected for light travel to the 
barycentre. A least-squares fit to all of these data yielded 
the ephemeris 



BMJD(TDB) = 53594.206270(1) +0.1582061029(3)^ 

for the white dwarf mid-egress. This ephemeris fits well 
with our 2004/2005 data (Fig. g}, but the historical data 
still show the ~100s deviat ion from a linear ephemeris dis- 
cussed by I Wolf et ail (|l993h . which they argued was caused 
by a third body in the system. An alternative explanation 
is that some of the historical timings are unreliable: the flux 
variations we observe in the egress feature (Figure [3]) sug- 
gest that some observations with small telescopes may have 
failed to detect the egress feature at all. This is particularly 
likely in the instances wh en a very long egress was reported 
(|Wood fc Crawf ord 1986), for which we see no evidence. 



5 DISCUSSION 

In this section we discuss the results of our model fits in more 
detail. The photometric method for parameter determina- 
tion in eclipsing CVs is dependent on precise measurement 
of the contact phases of the white dwarf and bright spot 
eclipses. These measurements are used to constrain the ge- 
ometry of the binary system, and derive the mass ratio and 



inclination. Our data affords the possibility of more precise 
determinations than previous authors, since we have for the 
first time identified the white dwarf ingress, and so can accu- 
rately determine the phase width of the white dwarf eclipse. 
We begin in Section [5.11 by discussing our determinations of 
the mass ratio and binary inclination, and the uncertainties 
on these measurements. We compare our results to deter- 
minations made both photometrically and spectroscopically 
by previous authors. In Section [5.21 we go on to investigate 
the individual components, and give masses, radii and ra- 
dial velocity semi-amplitudes for both the white dwarf and 
the donor star. In Section [5.31 we use these new parameter 
determinations to provide an updated estimate of the dis- 
tance to IP Peg. In Section 15.41 we use our determinations 
of the white dwarf colours to determine its temperature, by 
comparing our measurements to theoretical cooling models. 
This task is complicated by the variability we observe in 
the white dwarf egress feature. In Section 15.51 we detail a 
photoelectric absorption m odel which could account fo r this 
variability. As discussed in lTownslev fe Bildstenl (|2003h . the 
temperature of the white dwarf is driven by compressional 
heating of accreted matter, and so in Section f5. 61 we discuss 
the implications of our temperature measurement for the 
accretion history of this system. Finally, in Sections [5T7] and 
15.81 we examine the bright spot and accretion disc in more 
detail. In the case of the bright spot we determine its tem- 
perature and flux distribution, and for the accretion disc we 
examine changes in disc flux and radius over the course of 
our observations. 

5.1 Mass ratio and inclination 

The phase width of the white dwarf eclipse is an observable 
quantity that is intrinsically linked to two physical proper- 
ties: the mass ratio and the binary inclination. For a higher 
binary inclination the duration of the eclipse will be greater, 
thus to maintain the same phase width, as the inclination is 
increased, the size of the donor, and hence the mass ratio, 
must be decreased. There is th erefore a uni que relationship 
between these two properties ( Bailevl H979I ). This degener- 
acy can be broken since we have an additional geometric 
constraint due to the ingress and egress of the bright spot. 
The path of the accretion stream and hence the position of 
the bright spot is modified by the mass ratio. With this ad- 
ditional information we can determine the mass ratio and 
inclination in this system. 

Previous studies of this system did not identify the 
white dwarf ingress feature, and so there has been some un- 
certainty in the ph a se wid th A(j> of the white dwarf eclipse. 
IWood fc Crawford! (|l986l ) suggested that the white dwarf 
ingress was blended with the bright spot ingress, implying 
0.0863 < A(p < 0.0918. We find the phase width to be 
0.0935 ± 0.0003, implying either a higher mass ratio or a 
higher inclination than previously suggested, or some com- 
bination of the two. We find a much better fit for the high 
inclination case, due to the constraints implied by the bright 
spot. 

Our MCMC results for the mass ratio versus inclina- 
tion are plotted in Figure [5] The coloured contours and as- 
sociated points show the single parameter la errors and 
95% joint confidence contours. The other points plotted 
are various q and i determinations taken from the liter- 
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Figure 3. The region of the lightcurves around the white dwarf egress for sixteen of the eighteen eclipses. We label the eclipses according 
to the numbers given in the observation log (Table [TJ. We plot the u'-band data in the first column, the g'-band data in the second and 
the r'-band data in the third. For the x-axis timings we used the ephemeris given in Section I4.3I The dashed vertical lines identify the 
mean phase width of the white dwarf egress. 



ature (IWood fc Crawford! Il98rj ; iMarshl Il988l ; ISmakl |2002| ; 
iRibeiro et al.l I2007T ). The two black lines running across 
the plot show the uni que q/i relationship for ph ase widths 
of 0.0863 and 0.0935. IWood fc Crawford! i] 1986ft suggested 
0.0863 was the most likely value for the phase width: a num- 



ber of authors have assumed this value and hence underes- 
timated the inclination. 

If we now look at the values from our MCMC fits as 
plotted in Figure [5] it can be seen that the results from the 
different filters are not consistent at the 95% level. There 
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Table 2. White dwarf parameters obtained from MCMC fits to the individual lightcurves, as detailed in Section l4,2l We give the time of 
the white dwarf mid-egress tmid—egmss an d the white dwarf Ri and accretion disc radii (Ri and Rdisc) scaled by the binary separation 
a, as obtained from the g'-band fits. Rdisc = Rspot in these model fits. We list also the u 1 -, g'- and r'- band magnitudes of the white 
dwarf, determined from the step height of the white dwarf egress feature. 



Eclipse 


cycle 


^mid — egress 




ill /a 






Rspot j 





u' mag 






g' mag 






r' mag 


1 


-2200 


53246.145370(10) 





0037 ± 


0003 





262 ±0 


002 








18 


83 ±0 


02 








2 


-2194 


53247.094620(12) 





0043 ± 


0002 








19 


08 ±0 


08 


18 


53 ±0 


02 


19 


57 ±0 


06 


1 





53594.198840(8) 





0031 ± 


0008 





279 ±0 


002 


21 


29 ±0 


18 


18 


64 ±0 


01 


19 


21 ±0 


02 


5 


6 


53595.148038(5) 





0044 ± 


0011 





269 ±0 


001 


19 


50 ±0 


01 


18 


51 ±0 


01 


19 


19 ±0 


05 


6 


18 


53597.046564(5) 





0052 ± 


0008 





264 ±0 


002 


19 


51 ±0 


05 


18 


38 ±0 


01 


18 


64 ±0 


01 


7 


19 


53597.204771(6) 





0062 ± 


0008 





271 ±0 


002 


19 


14 ±0 


03 


18 


62 ±0 


01 


18 


89 ±0 


01 


8 


2 1 


53597.995814(5) 


o 


0060 ± 


0017 


o 


270 ± 


001 


19 


01 ±0 


03 


18 


36 ± 


01 


18 


61 ±0 


01 


9 


87 


53607.962749(5) 





0044 ± 


0005 





268 ±0 


003 


18 


73 ±0 


03 


18 


29±0 


01 


18 


54 ±0 


01 


10 


88 


53608.120980(8) 





0040 ± 


0006 








19 


15 ±0 


06 


18 


36 ±0 


01 


18 


84 ±0 


02 


11 


101 


53610.177655(4) 





0039 ± 


0007 





253 ±0 


001 


19 


63 ±0 


07 


18 


74 ±0 


01 


18 


23 ±0 


02 


12 


106 


53610.968692(5) 





0043 ± 


0005 





257 ±0 


001 


18 


91 ±0 


04 


18 


45 ±0 


01 


18 


87 ±0 


01 


13 


107 


53611.126905(5) 





0050 ± 


0008 





254 ±0 


001 


18 


47 ±0 


03 


18 


20 ±0 


01 


18 


41 ±0 


01 


14 


112 


53611.917924(2) 





0056 ± 


0010 





250 ±0 


002 


19 


15 ±0 


17 


18 


48 ±0 


02 


18 


94 ±0 


05 


15 


113 


53612.076134(4) 





0044 ± 


0005 





258 ±0 


001 


18 


64 ±0 


02 


18 


25 ±0 


01 


18 


65 ±0 


01 


16 


119 


53613.025377(7) 





0043 ± 


0006 





256 ±0 


001 


19 


08 ±0 


05 


18 


88 ±0 


01 


19 


23 ±0 


02 


17 


120 


53613.183571(7) 





0035 ± 


0007 





260 ±0 


001 


19 


57 ±0 


09 


19 


07 ±0 


01 


19 


56 ±0 


03 


18 


133 


53615.240256(10) 





0054 ± 


0010 





254 ±0 


002 


18 


86 ±0 


07 


18 


45 ±0 


01 


18 


65 ±0 


02 



-2200 



20 



40 



60 80 
Cycle no. 



100 



120 



140 



Figure 4. (O — C) values plotted against cycle number, using the 
eclipse timings listed in Tableland the linear ephemeris given in 
Section l4~3l We plot only our 2004/2005 datapoints, but used all 
available historical data in the determination of the ephemeris. 
Note the break in the x-axis of this plot between the 2004 and 
2005 data. 



are systematic uncertainties in our parameter determina- 
tions which dominate the uncertainties implied by the in- 
dividual fits. We believe this to be due to the uncertainty 
in the measurement of the phase widths. While the posi- 
tion of the white dwarf egress is well defined, the white 
dwarf ingress, and to some extent, the bright spot ingress 
and egress, are still affected by flickering, even in the phase- 
folded lightcurve. In order to quantify the uncertainty in our 
system parameters, it was necessary to estimate the range of 
variation in the phase width as caused by the flickering. Our 
phase- folded lightcurve is the average of the eleven complete 
lightcurves. If the position of the white dwarf ingress is af- 
fected by flickering, it follows that if we choose a subset of 
these eleven to make a new phase-folded lightcurve, then we 
will measure a different phase width. It is likely that in some 
lightcurves the flickering influences the apparent position of 
the ingress more than in others, and so if the phase width is 




0.3 0.4 0.5 0.6 

Mass ratio (q = M 2 /Mj) 

Figure 5. Mass ratio versus inclination, from our MCMC fits. 
We show the results of our fits to the phase-folded lightcurve, 
plotting the results in u' , g' and r' separately (blue, green and 
red, respectively). We plot points with the single parameter la 
errors and 95% joint confidence contours. The two solid lines run- 
ning diagonally across the plot indicate the q/i relationship for 
a phase width of 0.0863 (bottom) and 0.0935 (top). The inset 
shows an expanded view and the 68, 95 and 99% joint confidence 
contours. The grey points are measurements of q and i from the 
literature, and are label led a ccordin g to s ourc e as follows: W, 
[ Wood fc Crawford! lll986ri ; M. lMarshl (Il988h ; g. lSmakl J2002h ; R, 
iRibeiro et alJ l|2007l l. 
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measured in enough of these subsets, the uncertainty in the 
phase width can be determined, from which a more realistic 
uncertainty for the physical parameters can be calculated. 

We created a series of phase-folded lightcurves from dif- 
ferent combinations of the complete lightcurves. Each phase- 
folded lightcurve was made up of eleven lightcurves (like the 
original), but these el e ven were selected us i ng th e bootstrap 
method i|Efronl 1 19791 ; lEfron fc Tibshiranil 1 19931 ). For each 
of these a minimisation was performed using the MCMC 
method, as before, and the contact phases of the ingress 
and egress features were computed. We additionally calcu- 
lated the phase width of the white dwarf egress itself, using 
our previously determined q and i values, along with the 
weighted mean and variance of Ri/a calculated from the 
values listed in Table [2] We subsequently calculated the co- 
variance of the contact phases. These covariances are input 
into a Monte Carlo simulation through which we determined 
the uncertainty on the physical parameters of the system: 
q, i, the component masses M\ and M2, the binary separa- 
tion a and the radial velocity semi-amplitudes Kj and K^. 
One additional input into this calculation is an estimation 
of the extent to which the size of the white dwarf exceeds 
the Eggict on zero-temperature mass/radius relation (quoted 
in lVerbunt fc Rappaportlll988l ). For this the WD tempera- 
ture estimate fro m Section | 5.4 | is used with the white dwar f 
cool ing models o f Bergeron. Wese mael fc BeauchamrJ ()l995l ) 
and lHolberg fc Bergeror] (|2006l ). We find an oversize factor 
of ~1.04 to be appropriate. The results of this simulation are 
plotted in Figure [6] and all of these parameters are listed in 
Table [3] Unlike the earlier MCMC fits, the uncertainties on 
these determinations are more realistic values which account 
for the systematic errors due to the uncertainty in the phase 
widths, and so the determinations are consistent in all three 
bands. The (/-band determinations are of course the most 
precise, due to the higher signal-to-noise ratio in these data. 

We find the mass ratio in the (/-band data to be q — 
0.48 ± 0.01. This is consistent with, although more precise 
than, ma ny of the previous d eterminations: the photometric 
studies of lWood et al. 1119891 ) for example sugges ted a value 
of 0.49 , and lSmakl l|2002l ) reported q = 0.45+0.04. rWolf et all 
(|l993h suggested a value of 0.6 which is much higher than our 
finding, although previous photometric studies were limited 
by the difficulty in identifying the white dwarf ingress fea- 
ture in the li ghtcurves. However, the spectroscopic study of 
iMarshl l| 19881 ) also reported a higher value for q (0.58±0.06). 

The constraints o n the b inary inclination were reported 
by IWood fc Crawford! i|l986p to be 80.9 < i < 90 deg, and 
previous authors h ave su ggested values towards the low end 
of that range (e.g. ISmakll2002l . i = 81.8 ± 0.9 deg). Our g'- 
band value of 83.81 ± 0.45 deg is higher than most previous 
determinations. As we remarked earlier, a confounding fac- 
tor in previous studies has been the difficulty in identifying 
the white dwarf ingress feature, and the subsequent under- 
estimation of the phase width of the white dwarf's eclipse. 

5.2 The white dwarf and donor star parameters 

In Table [3] we list our determinations of the masses, radii 
and radial velocity semi-amplitudes of the two components, 
as determined from our Monte Carlo simulation. The white 
dwarf radius Ri is calculated from the weighted mean of the 
Ri/a values listed in Table [2j and the donor radius R2 is cal- 



cu lated using the Roche-lobe volume radius approximation 
of lEggletonl i|l983l ). 

If the (/-band results are used, the white dwarf mass is 
found to be 1.16 ± O.O2Af . This is larg e r tha n some pre- 
vious determinations, for example ISmakl (|2002l ) reported a 
value of 0.94 ± O.1OM , but is in very close agreement with 
IWatson et "aLI l|2003l ). who report a value of 1.16 - 1.18M . 
We find the donor mass to be 0.55 ± 0.0 2 M , which is 
greater than previous measure ments. [ Smakl J2002I) re ported 
a value of 0.42 ± 0.08 M , and IWatson et all (|2003h found 
M 2 to be O.5M . Our donor radius of O.47±O.O17? is con- 
sistent with previous authors, but it is low given our mass 
determination, suggesting a donor that is sligh tly undersized 
comp ared to an equivalent main sequence star (|Baraffe et alj 
119981 ). This is unusual for an accreting binary: donor stars 
tend to be oversized, due to either evolution or being out of 
thermal equilibrium as a result of mass loss. Our finding is 
consistent with a donor in thermal equilibrium. As we will 
show in Section [5.61 we find the accretion rate to be much 
lower than would be expected for this system, and this might 
explain why the donor is not significantly perturbed. If the 
donor is in thermal equilibrium then this suggests the mass 
transfer rate has been low for a very long period of time, 
in excess of the thermal timescale (~10 8 yr, for a star with 
these parameters). Alternatively, the mass transfer might 
have begun very recently. 



Recent work by iKniggel (|2006h has suggested a single, 
semi-empirical donor sequence for CVs with orbital peri- 
ods < 6h. This sequence suggests a donor mass of O.28M 
and radius of O.38i? for a P — 3.8h CV. These values are 
much lower than our measurements; for our donor mass of 
O.55M to be consistent with the sequence the binary period 
would need to be ~5.3h. IP Peg is not the only syst em that 
fits po orly to the semi-empirical donor sequence of IKniggel 
(2006); we see in figure 9 of that paper that there are a num- 
ber of other systems, such as RX J0944, which have a donor 
that is of an earlier spectral type than would be expected 
for their binary period. Again, the low accretion rate we will 
report in Section[5]6]may imply an unusual accretion history 
for this system, which could account for this discrepancy. 

Using our binary parameters, we calculated the ra- 
dial velocity semi-amplitudes K\ and K2, and find them 
to be 151 ± 3 and 317 ± 2 km/s respectively. Our value of 
K2 is larger than the long-established value (29 8 ± 8km/s; 
iMartin. Jones fc Smith! 19871 : iMartin et al1ll989l). wh ich was 



more recently supported by Beekman et al. (120001) . How - 



ever, the much better quality data of Watson et alj (|2003l ) 
give a value of 323.2 ± 3.5km/s. The consistency of our 
value derived purely from photometry with the best spec- 
troscopic measurement is a rare test of the consistenc y of 
the photometric method (see also iLittlefair et aill2008l and 
iTulloch. Rodriguez-Gil fc Dhillonl 2009 ) . The systematic un- 
certainty in spectroscopic determinations of Ki is much 
higher than in K2 , and so previous values of K\ in the litera- 
ture as determi ned from the emission lines vary significantl y 
(175±15 km/ s lMarshlll988l : 118±10 km/s, lHessmanlll989l ). 
More recently. ISmakl i|2002l ) used the three-body approxima- 
tion to find a value for K\ of 134±15 km/s. This is consistent 
with our finding of 151 ± 3 km/s. 
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Figure 6. The results of our Monte Carlo simulation, which we use to determine the parameter uncertainties listed in Table [3] We 
plot the mass ratio q versus the inclination i, the white dwarf mass Mi versus the donor mass M2 , and the white dwarf radial velocity 
semi-amplitude K\ versus the donor semi-amplitude K.%. The three columns separate out the results from the individual bands (left to 
right, u' , g' and r'). 



Table 3. Binary parameters for IP Peg. We list the phase width of the white dwarf eclipse, the mass ratio q, the inclination i and masses, 
radii and radial velocity semi-amplitudes of the two components. All of these parameters except the radii are determined from the Monte 
Carlo simulations described in Section I5.ll The white dwarf radius is derived from t he weighted mea n of the values determined in the 
individual MCMC fits, and the donor radius is calculated using the approximation of Egglcton (1983). 
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5.3 Distance to IP Peg 

Using our new parameters for the donor star, we can update 
the estimation of the distance to IP Peg. The surface b right- 
ness method of distance determination |Bailevl fl98ll ) uses 
the if-band magnitude to determine the distance to CVs. 
This method assumes all of the if-band flux originates in 
the donor star and is thus a lower limit, since the accretion 
disc will contribute an indeterminate amou nt of this flux (al- 
though this is small in the case of IP Peg) . ISzkodv fc Mated 
(| 1986T) measure d the K flux and estimated the distance to 
be 130 - 142pc. lBailevl |l98ll) assumed a constant donor sur- 
face brightness in his relation, but subsequent au thors have 
shown this quantity var ie s with the V — K colo ur (|Ramseverl 
1 1994 lBeuermannl(200rj) . iFroning et ail l| 19991 ) r eported an 
updat ed value of 121pc using this refine ment oflRamsevei 
1994) and the if-ba nd measurements o f ISzkodv fc Mateo 



1986). We applied the lBeuermarm] (|2006l ) update of Bailey's 



method, and fou nd the distance to IP Pe g to be 151pc. For 
this we used the lSzkodv fc Mated l|l986l ) -ff-band measure- 
ment, the donor r adius listed in T ableland the same V — K 
that was used bv lFroning et al] |l999), which is consistent 
with a donor spectral type of ~M4.8. The uncertainty in 
our distance determination is likely to be dominated by the 
uncertainty in V — K. We took ±1 to be an appropriate esti- 
mate of the uncertainty in the spectral type, and determined 
the su bsequent range inV — K using Table 3 of iBeuermannl 
(2006). We hence estimate the uncertainty in our distance 
determination of 151pc to be ±14pc. The best trigonometric 
estimate of the distance is 152+2gpc (Thorstensen, private 
communication), which is consistent with our result. 



5.4 White dwarf colours and temperature 

The main motivation for this work was a photometric deter- 
mination of the white dwarf temperature in this system. 
We see in Figure [3] that, while the length of the white 
dwarf egress appears constant, its flux varies significantly 
on timescales of days. Using the white dwarf magnitudes 
listed in Table [2] we calculate the u' — g' and g — r' colours 
of the white dwarf for each individual eclipse. We plot these 
data in the left panel of Figure [7J We plot also th e synthetic 
white dwarf tracks of lHolberg fc Bergeron! l|2006h for a white 
dwarf mass range of 0.3 to 1.2M0. 

This plot shows that as well as variation in all three 
bands, there are considerable variations in the white dwarf 
colours, particularly u — g' . There is a suggestion of an in- 
verse correlation between the u' — g and g' — r' colours. This 
is due to large variation in the u and r' fluxes and compar- 
atively little change in the g flux, so things become a little 
clearer when we plot u' — g' against the absolute (/-band 
magnitude M' g (Figure[7J right panel). We assume a distance 
of 151pc (Section |5.3[) . but plot the error bar corresponding 
to our estimate of the uncertainty in this determination. De- 
spite this uncertainty, we see in Figure [7] a clear trend in the 
data. The datapoints with the lowest u — g' colours inter- 
sect with the synthetic tracks, and there are many other 
points which are much redder in colour. Not only are many 
of these points too red to correspond to a white dwarf of 
mass ~1.1M0, they are too red to correspond to a white 
dwarf of any mass. 

We consider two possibilities for the variation in the 



white dwarf colours. Firstly, it is possible that there is 
some variable contamination from an additional tempera- 
ture component, such as a band of hot, a ccreted material 
aroun d the equator of the white dwarf (|Piro fc Bildstenl 
l2004h . If this was the case one expects the trend in the 
egress datapoints to be approximately in line with the 
direction of increas i ng wh ite dwarf temperature in the 
iHolberg fc Bergeron! ()2006l ) tracks. However, what we ac- 
tually see is a trend that is roughly perpendicular to the di- 
rection of increasing white dwarf temperature. The inverse 
correlation between the u'—g' and g'—r' colours is due to the 
fact that there is comparatively little change in the g' flux 
compared to the (correlated) change in the u' and r' fluxes. 
It is difficult to reconcile this with a model that includes an 
additional temperature component. It is more likely that the 
variation in the white dwarf colours is due to photoelectric 
absorption, as a result of obscuration of the white dwarf by 
material from the accretion disc. This model can approxi- 
mately explain the colour variations, and we examine this in 
more detail in Section [53] An occulting region of disc mate- 
rial could explain past r eports of very extended w hite dwarf 
egresses m this system (|Wood fc Crawford! Il986h . Alterna- 
tively these might have been failures to detect the variable 
egress at all. Given the short period of time over which we 
observed this system it is difficult to quantify the extent of 
this phenomenon. 

We now attempt to estimate the white dwarf tem- 
perature in this system via comparison with the synthetic 
tracks. Figure [7] suggests a temperature of between 10,000 
to 15, 000K; it is difficult to be more precise than this due 
to the flux variations. Even so, this range is remarkably 
cool. Previous studies have suggested that a white dwarf 
in a CV with a period > 3h would h ave a typical temper- 
ature of between 25, 000 - 35,000K (jWinter fc Sionl 120031 ; 
iTownslev fc Bildstenl lioolh . Such a temperature is clearly 
inconsistent with our IP Peg data: this would require the 
white dwarf to be more than a magnitude more luminous in 
g' and more than 1.5 magnitudes brighter in u . In order to 
plot the absolute magnitudes used in the right panel of Fig- 
ure [7] we assumed a distance to IP Peg and the uncertainty 
in this measurement could allow a higher temperature, but 
judging from the parallax measurement this uncertainty is 
unlikely to be sufficient to reconcile a 25, 000K white dwarf 
with the observations. In addition, the v! — g' / g' — r' plot 
is independent of distance and also implies a temperature 
range of 10,000 to 15, 000K for the white dwarf. 

In making this temperature determination we have as- 
sumed that our points with the bluest u' — g colour are 
indicative of a white dwarf that is largely unobscured. How- 
ever, we also investigate the possibility that the accretion 
disc is opaque, blocking the bottom half of the white dwarf 
from view and so producing a situation where the obscura- 
tion of the white dwarf is at least 50% at all times. With this 
modification we find the white dwarf temperature implied 
by our model fits is at most ~21,000K. This is higher than 
our previous findings, but still significantly smaller than the 
predicted 25, 000 - 35, 000K. Additionally, we will show in 
Section 15.51 that the colour trends in the datapoints can be 
explained by an absorption model if we assume a 15, 000K 
white dwarf, but not if we assume a temperature of 20, 000K 
or higher; hence we prefer the fits which imply a cooler white 
dwarf. 
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Figure 7. Using the eclipse determinations listed in Table [2] we plot datapoints showing the colours and magnitude of the white dwarf 
for each individual eclipse. In the left panel we plot the u' — g 1 and g' — r' colours for the white dwarf, and in the right panel we show the 
u' — g' colour versus the absolute magnitude M' g . For this we use our distance to IP Peg of 151pc (Section 15. 31 . but we indicate to the 
right of the plot our estimate of th e uncertainty in the d i stanc e modulus. The red lines that run from the top left area of each panel are 
the synthetic white dwarf tracks of lHolbere &i Bergeron! J2006j). In the left panel we plot tracks for a white dwarf mass of 0.3 to 1.2Mq, 
plotted every O.IMq, but for clarity in the right panel we just show the tracks for a mass of 1.0, 1.1 and 1.2Mq. Perpendicular to these 
tracks we plot isotherms (blue) every 5, 000K. The green lines are representative of the results of our opacity model, discussed in more 
detail in Section 15.51 Line (a) assumes a white dwarf with a temperature of 15,000K and a mass of I.IMq, and absorbing material with 
a temperature of 8,000K and Nh < 10 23 cm — 2 . Line (b) uses the same values as line (a), except for the absorber temperature which 
is assumed to be 14,000K. Line (c) uses the same values as line (a), except for the white dwarf temperature which is assumed to be 
20,000K. 



5.5 An absorption model 

We propose that the variation in the white dwarf colours is 
due to occultation of the white dwarf by disc material, caus- 
ing a significant and varying photoelectric absorption of the 
white dwarf flux. To test this we computed an opacity model 
and compared it to the data. Our intention here is not to 
make a precise and quantitative description of the absorp- 
tion in this system, but to determine if a varying opacity is 
an adequate description for the flux variations we observe. 

We calculated the opacity for a range of absorber tem- 
peratures and neutral hydrogen column densities, allowing 
for hydrogen bound-free, free-free, bound-bound and H~ 
opacities. A white dwarf spectrum was used to determine the 
source flux and we assumed LTE, an electron density of 10 14 
cm* 3 , and we approximated the line broadening by assum- 
ing it to be Gaussian with an RMS of 50 km/s. We convolved 
the resultant spectra with the SDSS filter profiles and hence 
calculated the change in the u' , g' and r 1 magnitudes for 
each absorber temperature and column density. These cal- 
culations give us the change in the source flux; we chose the 
unabsorbed colours to be those of a 15, 000K, I.IMq white 
dwarf, which is appropriate given our earlier measurements. 



The results of this model calculation are plotted in Figure [7J 
as tracks (a) and (6). These two tracks denote a constant ab- 
sorber temperature of 8, 000 and 14, 000K respectively. The 
hydrogen column density increases along these tracks. If we 
compare the tracks to the white dwarf measurements we see 
that they lie approximately in line with the trend we observe 
in the datapoints, suggesting that absorption is an adequate 
explanation for the variations in white dwarf flux. There are 
some discrepancies; in particular a number of datapoints at 
low temperatures are not consistent with the tracks we plot, 
but we attribute this to the simplicity of our model. 

It was noted in Section 15.41 that it is possible that in 
all of our observations we see only a small fraction of the 
white dwarf flux, which would make our determination of 
the white dwarf temperature a significant underestimation. 
To test this we plot a third model track in Figure [7] la- 
belled (c). This track uses the same parameters as track (a) 
with the exception of the white dwarf temperature, which 
is assumed to be 20, 000K. We see this track is a poor fit 
to our observations, and so a white dwarf temperature of 
20, 000K or greater is not supported by our opacity model. 
In addition, if the temperature of the white dwarf is greater 
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than ~17, OOOK, we find that there is no adjustment in the 
temperature of the absorbing material that can lead to the 
redwards movement in g — r that would be required to fit 
the data. 

Based on the neutral hydrogen column densities we use 
in our model, we calculate an Nh of up to ~10 23 cm -2 is nec- 
essary to account for the v ariation we observe in the white 
dwarf flux. For comparison iHorne et alj (| 1994T ) reported ab- 
sorption of the white dwarf emission by disc material in the 
dwarf nova OY Carinae, but calculated an Nh value that is 
an order of magnitude lower than we find for IP Peg. The 
higher Nh value we find for IP Peg may be due to the very 
high inclination of this system. If this value is correct, it will 
be extremely difficult to detect the white dwarf in IP Peg 
at ultraviolet wavelengths. The ultraviolet co ntinuum out- 
of-ec l ipse flux has been measured i n IP Peg (|Hoard et al.l 
ll997l ; ISaito. Baptista fc Hornell2005h , and is consistent with 
a cool white dwarf and/or heavy absorption of the white 
dwarf flux. 

5.6 The long-term accretion rate 

iTownslev fc Bildstenl (|2003h noted that the temperature of 
accreting white dwarfs is a good tracer of the medium term 
(10 3 -10 5 yr) accretion rate. Our value for the IP Peg white 
dwarf implies a mean accretion rate of < 5 x 1O _11 M0 
yr -1 . The expected rate for a dwarf nova of this period is 
~2 x 1O -9 M0 yr -1 : more than 40 times greater. The value 
we find for the mean accretion rate is much closer to the in- 
stantaneous rate, which we estimate to be ~1.8 x 1O _1O M0 
yr -1 by following the calculation of iMarshl (|l988T ) and sub- 
stituting in our new distance determination. 

Our determination of the mean accretion rate implies a 
mass transfer timescale of ~10 10 yr, much greater than the 
thermal timescale for the donor, which is ~10 8 yr. We would 
therefore expect the donor in IP Peg to be unperturbed by 
the mass transfer, in contrast to the donors in many other 
CVs which are observed to be oversized when compared to 
isolated main sequence stars of equivalent mass. As we noted 
in Section 15.21 we do find a mass and radius for the donor 
that are consistent with those of a main sequence star in 
thermal equilibrium. 

The low temperature of this system suggests there could 
be a selec tion bias in the sampl e of lon g-period systems pre- 
sented in lTownslev fc Bildstenl (|2003l ). and further detailed 
studies of eclipsing systems are necessary in order to de- 
termine the extent of this bias. At this stage, for example, 
it is possible that mass transfer in IP Peg has begun very 
recently and so the temperature in this system is a poor 
indicator of the long term rate. One p ossible mechanism i s 
the 'hibernation scenario' proposed by IShara et al.l <|l986h . 
in which, following a classical nova eruption and a white 
dwarf cool-down period, the binary is detached for 10 3 - 
10 6 yr. |Sh ara et al. I i|l986l ) argue that an old nova emerging 
from this period of detachment (due to angular momentum 
loss bringing the donor back into contact with its Roche 
lobe) might first appear as a dwarf nova. 

Alternatively, we can take our low value for the mass 
transfer rate over the past 10 3 -10 4 yr in this system to be 
correct. This is surprising, given the contention that long- 
period dwarf novae have a high mean M despite their ob- 
served instantaneous values. The low temperature of IP Peg 



would suggest it has sustained an accretion rate well be- 
low the expected mean rate for more than 1000 years. If all 
'steady-state' CVs go through low accretion rate dwarf nova 
phases, then this result implies the length of these phases 
can be very long indeed. 

The third possibility is that there are problems with 
the classical picture of CV evolution as driven by an- 
gular moinerrtu^Jos^jrhe orbital period distribution of 
CVs l|Spruit fc R ittcr 1983) has traditionall y been explained 
by th e disrupted magnetic braking model (|Robinson et al.l 
Il98lf ). whereby angular momentum loss above the 2 - 3 h 
'period gap' is due to a combination of gravitational radia- 
tion and magnetic stellar wind braking, and below the gap is 
due to gravitational radiation only. It follows that angular 
momentum loss is higher at longer periods, and since an- 
gular momentum loss is proportional to M, M is expected 
to be high at long periods. The existence of low M dwarf 
novae like IP Peg at long periods conflicts wi th this picture. 
Some authors suggest alternatives, such as IWillems et al.l 
j2007l ). who consider the effect of circumbinary discs on 
C V evolution. The the oretical orbital period distribution 
of IWillems etafl (|2007h is a reasonable reproduction of the 
observed distribution, and they find that systems with un- 
evolved donors need not detach and evolve below the period 
gap as in the disrupted magnetic braking model. In their 
model different CVs can have very different long-term Ms, 
and some systems can exist as dwarf novae throughout their 
lifetime. This would explain the M we measure in IP Peg, 
although study of further systems is necessary to distinguish 
between this and the other possibilities. 

5.7 The bright spot 

It is apparent from the lightcurves (Figure [1} that most of 
the flickering originates in the bright spot. The bright spot 
is modelled as a linear region which approximately follows 
the outer circumference of the disc, with a luminosity which 
decreases approximately exponentially from the point of im- 
pact of the accretion stream (see Appendix The geome- 
try is such that it is the impact region which is eclipsed first, 
and is first to emerge from eclipse. In Figure [8] we plot the 
residuals from the model fits to the individual lightcurves, 
and we indicate on this figure the beginning and end points 
of the bright spot ingress and egress features. It is appar- 
ent in this figure that the majority of the variation due to 
flickering is localised in the impact region, since we see the 
residuals very quickly flatten off following the beginning of 
the ingress, but return immediately following the beginning 
of the egress. 

In Figure[9]we plot the colour variation along the bright 
spot. For this we used the bright spot parameters listed in 
Table lATI with Equation lA3l to calculate analytically the flux 
in each filter at each point along the length of the spot. We 
plot the bright spot colour from the point of impact of the 
gas stream on the accretion disc (labelled 'A' in Figure [9]l 
to a point 10 scale lengths (I; see Appendix I A4|l away from 
the impact region (labelled 'B'). We plot also a line showing 
the colours of blackbodies with temperatures in the ~6, 000 
- 17, OOOK range. At the point of impact the bright spot 
is close in colour to a 13, OOOK blackbody, and as we move 
away from the impact region the colours of the bright spot 
decrease approximately in the direction of decreasing black- 
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Figure 8. The residuals for the model fits to the individual </- 
band lightcurves. The four vertical lines show (from left to right) 
the beginning and end points of the bright spot ingress, followed 
by the begining and end points of the bright spot egress. A phase 
of zero corresponds to the mid-point of the white dwarf egress. 
We omit eclipse 2 because it includes no data for the bright spot, 
and eclipse 3 because the flare event near the white dwarf egress 
precludes a precise fit. 



body temperature. The spectroscopic study of lMarshl (Il98cf ) 
also found the closest blackbody fit to the bright spot was 
for a temperature of 13, 000K. The shape of the bright spot 
line is not very important since it is very model dependent 
and as we note in Appendix I A4I our bright spot model does 
not describe the IP Peg bright spot perfectly. However, the 
fact that the general trend of the line is in the direction 
of decreasing blackbody temperature is clear evidence for 
cooling along the line of bright spot elements. The variable 
temperature of the bright spot in IP Peg will cause it to vary 
according to the wavelength of observation. For instance 
near- infrared data suggest a cooler 10, 000K (|Froning et al.l 
1999), and hotter parts than we see may dominate in the 
far-ultraviolet. 

The points 'A' and 'B' are also labelled in the schematic 
of the primary Roche lobe. An additional point to note from 
the schematic is that the angle of the bright spot line of ele- 
ments is such that it runs approximately along the circum- 
ference of the accretion disc. This behaviour is not enforced 
in our model, but seems physically appropriate: the incom- 
ing material from the gas stream is thermalised at point 'A' 
when it strikes the dense accretion disc, and as it cools, it 
propagates around the disc along a line of approximately 
constant radius and disc density. 



5.8 The accretion disc 

5.8.1 Accretion disc radius 

Data collected by the American Association of Variable Star 
Observer^ shows that IP Peg was observed to go into out- 
burst approximately 35 days prior to the beginning of our 
August 2005 observations. In Figure[l0]we plot the accretion 
disc radius scaled by the binary separation against the MJD 
time since the detection of the outburst, using the (/-band 
model fits. We see there is a decline in disc radius over the 
course of the ~20 days covered by our observations. During 
outburst the accretion disc radius increases as angular mo- 
mentum is transported outwards due to the disc viscosity, 
and following the outburst the disc radius begins to shrink 
as low angular momentum material is added to the disc via 
the gas stream. It is this decline we observe in Figure 1101 
We compared these val ues with prev i ous d eter minations of 
the di sc radius given bvlWood et al.l ljl989T ) and I Wolf et al.l 
l|l993h . In Figure 6 of I Wolf et all 4 1993 ) disc radius is plotted 
as a function of time since outburst, and our measurements 
are consistent with the rate of decline shown in their figure. 



5.8.2 Accretion disc flux 

The current and well-established paradigm for outbursts 
in dwarf novae is the therm al- viscous accretion disc in- 
stability model l|Osakil Q.974). This model proposes that 
the outbursts are due to the bi-stable nature of accre- 
tion discs when Tdisc ~10, 000K, at which point hydro- 
gen in the disc changes from a n eutral to an ionised state 
ijMever fc Mever-Hof mcistcr 1981). During quiescence, mat- 
ter is stored in the disc and at a critical amount the in- 
stability sets in and a significant fraction of this mass 
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Figure 9. Left: The solid line (black) shows the colour variation along the bright spot, from the point of impact of the gas stream 
(labelled 'A'), to a point 10 scale lengths away from this impact region (labelled 'B'). The dashed line (blue) shows for comparison the 
blackbody colours for a temperature ranging from ~6,000 - 17, 000K. Right: Schematic of the primary Roche lobe. The dashed line 
shows the extent of the Roche lobe, the dotted line shows the accretion disc radius and the central point is the white dwarf (all to scale). 
The line originating at the LI point is the gas stream, and the bright spot is plotted from the impact region 'A' to 'B' as a series of 
circles, with the area of each circle representing the surface brightness at that point. 
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and the decrease in flux from this region counterbalances 
the increase in the effective tem perature of the outer disc 
i|Truss. Wvnn fc Wheatlevll2004 ). 
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Figure 10. Variations in accretion disc radius over the course of 
the August/September 2005 observations. We plot the accretion 
disc radius Rdisc scaled by the binary separation against the MJD 
time since the last outburst. 



is very rapidly accreted onto the central object. There 
are a number of detailed one-dimen sional calculations of 
the disc instability mode l (see, e.g.. lHameurv et al.lll998l ; 
iBuat-Menard et al.l l200ll ). and among their predictions is 
a significant increase in the disc luminosity during quies- 
cence, since as matter is stored in the disc the surface den- 
sity increases, prompting an increase in the effective tem- 
perature. However this is not observed: long term studies 
of dwarf novae show th e visible flux is approxim ately con- 
stant during quiescence |Cannizzo fc M attci 1992). Much of 
the existing data has been collected by amateurs and is of 
variable quality, but the models predict an increase of ~1 
-3 visual m agnitudes, w hich should be easily detectable in 
these data (Smak 2000). It has been proposed that these 
observations could be explained by a hot, high-viscosity in- 
ner region of the disc which slowly cools during quiescence, 



Eclipsing systems such as IP Peg allow for a more pre- 
cise measurement of disc flux variations, since during the 
phase range over which both the white dwarf and the bright 
spot are eclipsed, the only contributors to the observed flux 
are the donor star and the accretion disc. The donor star 
component should be essentially constant, and so by mea- 
suring the average flux at a constant phase for every eclipse, 
we can determine the luminosity evolution of the accretion 
disc. We chose to measure the flux at the phase of the white 
dwarf egress, since at this point we observe a complete half 
of the disc and so our results should be unbiased by tempo- 
ral radial variations. We used our model fits to subtract the 
white dwarf contribution at this point. We plot the average 
flux over a 10s period at mid-egress in Figure 1111 and see 
an apparent decrease in disc flux over the period covered, 
in contrast to the increase which is predicted by the disc 
instability model. Our findings are c onsistent with the 'mir - 
ror eclipse' phenomenon reported bv lLittlefair et all l|200ll ). 
which they interpreted as eclipses of the donor star by the 
optically thin outer edge of the accretion disc. This conflicts 
with the disc instability model, which requires an optically 
thick disc for physically viable values of the disc viscosity 
parameter. Note also that the decrease in brightness we ob- 
serve (from ~0.4 to O.lmJy) cannot be explained purely in 
terms of the disc radius variations we discussed in Section 
15.8.11 . since we see in Figure [10] that the change in disc area 
over the course of our observations is only ~25%. Our results 
demonstrate a genuine dimming of the disc over the course 
of our observations. We also do not ob serve any ev i dence for 
the hot inner disc region proposed by [Truss et all (I2004T I. 
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Figure 11. Variations in accretion disc flux over the course of 
the August/September 2005 observations. We plot the (/-band 
accretion disc flux at the point of the white dwarf egress against 
the MJD time since the last outburst. 
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6 CONCLUSIONS 

In 2004 and 2005 we took eighteen separate observations 
of the eclipse in the 3.8h period dwarf nova IP Peg. We 
used the triple-band camera ULTRACAM mounted on the 
William Herschel Telescope, and observed simultaneously 
in the Sloan u'-, g'- and r'-bands. By phase-folding these 
data we identify for the first time the white dwarf ingress 
in this system. The phase width of the white dwarf eclipse 
is 0.0935 ± 0.0003, significantly higher than the standard 
value used by previous authors. This new value implies a 
higher orbital inclination or a higher mass ratio or some 
combination of the two. Our model fits strongly favour a 
high orbital inclination, which we find to be 83.81T0.45 deg. 
Our determinations of the masses of the two components are 
higher than those reported by previous authors. 

We find the brightness of the white dwarf egress feature 
to be highly variable in flux and colour. We suggest this is 
due to photoelectric absorption as a result of obscuration of 
the white dwarf by material from the accretion disc. We find 
the duration of the egress feature to be consistent in length, 
and so suggest that past reports of very extended egresses 
were due to failures to detect the egress at all, as a result of 
this varying obscuration. 

It is difficult to determine an accurate temperature for 
the white dwarf due to the obscuration, but by comparing 
our data with synthetic white dwarf models we find the most 
likely temperature to be in the range of 10, 000 - 15, 000K. 
This is v ery low for a dwarf nova abo ve the period gap and, 
following iTownslev fe Bildstenl (|2003T ). implies a mean accre- 
tion rate over the past 1, 000 - 10, 000 years of < 5x 1O _11 M 
yr -1 . This is more than 40 times lower than the expected 
rate. Additionally we find a mass and radius for the donor 
that are consistent with those of a main sequence star in 
thermal equilibrium, which also suggests a very low accre- 
tion rate. Unless the mass transfer in IP Peg has begun very 
recently, these findings imply either that CVs can sustain 
accretion rates well below the expected rate for very long 
periods of time, or that the classical picture whereby the 
long-term accretion rate scales with period is not correct. 



REFERENCES 

Bailey J., 1979, MNRAS, 189, 41P 
Bailey J., 1981, MNRAS, 197, 31 

Baraffe I., Chabrier G, Allard F., Hauschildt P. H., 1998, 
A&A, 337, 403 

Beekman G, Somers M., Naylor T., Hellier C, 2000, MN- 
RAS, 318, 9 

Bergeron P., Wesemael F., Beauchamp A., 1995, PASP, 

107, 1047 
Beuermann K., 2006, A&A, 460, 783 

Buat-Menard V., Hameury J.-M., Lasota J.-P., 2001, A&A, 
366, 612 

Cannizzo J. K., Mattei J. A., 1992, ApJ, 401, 642 

Dhillon V. S., Marsh T. R., Stevenson M. J., Atkinson 
D. C, Kerry P., Peacocke P. T., Vick A. J. A., Beard 
S. M., Ives D. J., Lunney D. W., McLay S. A., Tierney 
C. J., Kelly J., Littlefair S. P., Nicholson R., Pashley R., 
Harlaftis E. T., O'Brien K., 2007, MNRAS, 378, 825 

Efron B., 1979, Annals of Statistics, 7, 1 

Efron B., Tibshirani R. J., 1993, An Introduction to the 
Bootstrap. Chapman & Hall, New York 

Eggleton P. P., 1983, ApJ, 268, 368 

Feline W. J., Dhillon V. S., Marsh T. R., Brinkworth C. S., 

2004, MNRAS, 355, 1 
Feline W. J., Dhillon V. S., Marsh T. R., Watson C. A., 

Littlefair S. P., 2005, MNRAS, 364, 1158 
Froning C. S., Robinson E. L., Welsh W. F., Wood J. H., 

1999, ApJ, 523, 399 
Gansicke B. T., Beuermann K., de Martino D., 1995, A&A, 

303, 127 

Goranskii V. P., Lyutyi V. M., Shugarov S. Y., 1985, Soviet 

Astronomy Letters, 11, 293 
Goranskii V. P., Shugarov S. Y., Orlowsky E. I., Rahimov 

V. Y., 1985, Information Bulletin on Variable Stars, 2653, 

1 

Hameury J.-M., Menou K., Dubus G, Lasota J.-P., Hure 

J.-M., 1998, MNRAS, 298, 1048 
Hessman F. V., 1989, AJ, 98, 675 



© 0000 RAS, MNRAS 000, 000-000 



Physical properties of IP Pegasi 17 



Hoard D. W., Baptista R., Eracleous M., Home K., Mis- 
selt K. A., Shafter A. W., Szkody P., Wood J. H., 1997, 
MNRAS, 288, 691 

Holberg J. B., Bergeron P., 2006, AJ, 132, 1221 

Home K., Marsh T. R., Cheng F. H., Hubeny I., Lanz T., 
1994, ApJ, 426, 294 

Knigge C, 2006, MNRAS, 373, 484 

Littlefair S. P., Dhillon V. S., Marsh T. R., Gansicke B. T., 

2006, MNRAS, 371, 1435 
Littlefair S. P., Dhillon V. S., Marsh T. R., Giinsicke B. T., 

Southworth J., Baraffe I., Watson C. A., Copperwheat C, 

2008, MNRAS, 388, 1582 

Littlefair S. P., Dhillon V. S., Marsh T. R., Harlaftis E. T., 

2001, MNRAS, 327, 475 
Lubow S. H., Shu F. H., 1975, ApJ, 198, 383 
Marsh T. R., 1988, MNRAS, 231, 1117 
Martin J. S., Jones D. H. P., Friend M. T., Smith R. C, 

1989, MNRAS, 240, 519 
Martin J. S., Jones D. H. P., Smith R. C, 1987, MNRAS, 

224, 1031 

Meyer F., Meyer-Hofmeister E., 1981, A&A, 104, L10+ 

Osaki Y., 1974, PASJ, 26, 429 

Piro A. L., Bildsten L., 2004, ApJ, 610, 977 

Press W. H., 2002, Numerical recipes in C++ : the art of 
scientific computing. Cambridge University Press 

Pyrzas S., Gansicke B. T., Marsh T. R., Aungwerojwit A., 
Rebassa-Mansergas A., Rodriguez-Gil P., Southworth J., 
Schreiber M. R., Nebot Gomez-Moran A., Koester D., 

2009, MNRAS, 394, 978 
Ramseyer T. F., 1994, ApJ, 425, 243 

Ribeiro T., Baptista R., Harlaftis E. T., Dhillon V. S., Rut- 
ten R. G. M., 2007, A&A, 474, 213 

Ritter H., Kolb U., 2003, A&A, 404, 301 

Robinson E. L., Barker E. S., Cochran A. L., Cochran 
W. D., Nather R. E., 1981, ApJ, 251, 611 

Saito R. K., Baptista R., Home K., 2005, A&A, 433, 1085 

Shafter A. W., 1992, ApJ, 394, 268 

Shara M. M., Livio M., Moffat A. F. J., Orio M., 1986, 

ApJ, 311, 163 
Smak J., 2000, New Astronomy Review, 44, 171 
Smak J. I., 2002, Acta Astronomica, 52, 189 
Smith J. A., Tucker D. L., Kent S., Richmond M. W., 
Fukugita M., Ichikawa T., Ichikawa S.-i., Jorgensen A. M., 
Uomoto A., Gunn J. E., Hamabe M., Watanabe M., Tolea 
A., Henden A., Annis J., 2002, AJ, 123, 2121 
Southworth J., Hickman R. D. G., Marsh T. R., Rebassa- 
Mansergas A., Gansicke B. T., Copperwheat C. M., 
Rodriguez-Gil P., 2009, MNRAS, submitted, 
Spruit H. C, Ritter H., 1983, A&A, 124, 267 
Szkody P., Mateo M., 1986, AJ, 92, 483 
Townsley D. M., Bildsten L., 2003, ApJL, 596, L227 
Townsley D. M., Gansicke B. T., 2009, ApJ, 693, 1007 
Truss M. R., Wynn G. A., Wheatley P. J., 2004, MNRAS, 
347, 569 

Tulloch S. M., Rodriguez-Gil P., Dhillon V. S., 2009, MN- 
RAS, 397, L82 

Verbunt F., Rappaport S., 1988, ApJ, 332, 193 

Warner B., 1995, Cataclysmic variable stars. Cambridge 
Astrophysics Series, Cambridge, New York: Cambridge 
University Press 

Watson C. A., Dhillon V. S., Rutten R. G. M., Schwope 
A. D., 2003, MNRAS, 341, 129 



Willems B., Taam R. E., Kolb U., Dubus G., Sandquist 

E. L., 2007, ApJ, 657, 465 
Winter L., Sion E. M., 2003, ApJ, 582, 352 
Wolf S., Mantel K. H., Home K., Barwig H., Schoembs R., 

Baernbantner O., 1993, A&A, 273, 160 
Wood J., Crawford C. S., 1986, MNRAS, 222, 645 
Wood J. H., Marsh T. R., Robinson E. L., Stiening R. F., 

Home K., Stover R. J., Schoembs R., Allen S. L., Bond 

H. E., Jones D. H. P., Grauer A. D., Ciardullo R., 1989, 

MNRAS, 239, 809 



APPENDIX A: THE LIGHT-CURVE 
MODELLING CODE 

In this section we provide a concise description of the code, 
LCURVE, that we wrote in order to model th e light-curves 
of this pap er; the same code has bee n used in lPvrzas et al.l 
l|2009h and ISouthworth et all (|2009l 'l. The aim of the code 
was to provide a flexible framework for fitting light-curves 
characteristic of eclipsing dwarf novae and detached white 
dwarf / M dwarf binary stars. LCU RVE was developed as 
a generalisation of the code used in iHorne et all l|l994h to 
model the dwarf nova OY Car and subsequently employed 
in the an alysis of ULTRACAM photometry of several other 
systems (|Feline et al.ll2004 120051 ; iLittlefair et~aLll2006l ') . 

The binary is defined by up to four components: white 
dwarf, secondary star, accretion disc and bright-spot. In the 
case of detached systems, the latter two are disabled. Each 
component is described by a set of small, fiat elements, each 
of specified area, position, orientation and brightness. The 
light-curves are calculated by summing the contributions 
from all faces which are oriented towards the observer, allow- 
ing for eclipses by the two stars. Additionally, it is possible to 
make the accretion disc opaque or transparent. The contri- 
bution of a face depends upon how it is oriented with respect 
to the line of sight. Denoting the angle between the perpen- 
dicular to a face and the vector pointing towards Earth as 9, 
there is no contribution if 6 > 7r/2 radians, or equivalently 
if jj, = cos# < 0. Otherwise the contribution is typically 
specified using a limb-darkened Lambert's law behaviour: 

I <x n(l-e + en), (Al) 

where e is a linear limb-darkening coefficient (the code also 
allows for quadratic limb-darkening). 

The effect of finite length exposures is allowed for by 
calculating multiple points covering an exposure and using 
trapezoidally-weighted averaging. Thus 30 second-long ex- 
posures of a white dwarf eclipse can in the model be split 
into 7 points separated by 5 seconds for example. 

In lHorne et al.l (1 19941 ) , the overall scaling of the various 
components was handled using singular- value decomposition 
(SVD), allowable if all components contribute linearly to the 
light-curve. The advantage of this approach is speed: during 
the fitting process, several potential free-parameters are im- 
mediately removed at an early stage, allowing more time for 
optimisation of the remaining model parameters. However, 
it is not always a valid approach, in particular if "reflec- 
tion" (more properly reprocessing) is significant, then one 
cannot assume a simple linear combination, since doubling 
the contribution of the white dwarf, say, may also imply an 
altered contribution from the heated secondary star. Thus 



© 0000 RAS, MNRAS 000, 000-000 



18 CM. Copper-wheat et al. 



in LCURVE one can switch off the SVD scaling, falling back 
on the physical model to handle the effects of reflection. 

We experimented with several forms of minimis ation: 
the simplex method implemented via "amoeba" l)Pressl 
2002), Powell's method and the Levenburg-Marquardt 
method. The simplex method has the advantage of robust- 
ness, while Levenburg-Marquardt can deliver uncertainty es- 
timates once it has reached a minimum. However we found 
neither to be particularly secure in reaching the minimum 
\ 2 for the rather complex fits that we undertook in this pa- 
per. Early on in the fitting we realised on several occasions 
that we were trapped in local-minima, but even when we 
were not, the routines would often report that the fits were 
converged when subsequent work revealed this not to be the 
case. We believe that this was largely caused by degeneracies 
in the model which are very hard to remove entirely. There is 
for instance a fundamental degeneracy between mass ratio q 
and orbital inclination i that cannot be avoided. Ultimately, 
the method that helped above all in obtaining the final fits, 
and which, after extensive use, gave us confidence that we 
did achieve global minima, was the Markov Chain Monte 
Carlo (MCMC) method. Our standard technique evolved to 
become (i) carry out an initial minimisation with the sim- 
plex method, and (ii) then run an MCMC chain until it 
reaches an essentially flat run of \ 2 with fit number, con- 
sistent with the statistical fluctuations characteristic of the 
MCMC method. For our main fits we performed the same 
process from 10 independent starting models to check that 
we obtained consistent results. We typically required chains 
of 10,000 fits or more to reach the desired quasi-static state. 
We strongly recommend this technique in all similar work 
to guard against premature claims of "convergence". 

Apart from the general physical parameters mentioned 
above and others associated with the individual components, 
the other parameters needed to define the model, were com- 
putational, such as the number of elements needed to define 
the various components. Our strategy with these was sim- 
ply to make them large enough that they had no significant 
effect upon the best-fit parameters. We finish by describ- 
ing the model used for the individual components in more 
detail. 



Al The white dwarf 

We modelled the white dwarf as a linearly limb-darkened 
sphere. Although we also enabled the possibility of Roche- 
distortion, we never used it in the fits presented here. The 
eclipse of the three other components by the white dwarf 
was included by default, although its effect is tiny in the 
case of IP Peg since the white dwarf is so small. 



A2 The secondary star 

The secondary star was modelled with full accounting for 
distortion in Roche geometry. Again linear limb-darkening 
was assumed, and heating by flux from the white dwarf was 
included along with gravity darkening. None of this is es- 
pecially significant for IP Peg as the secondary star con- 
tributes relatively little. The important feature is of course 
the eclipse, again calculated accounting for Roche geometry. 



A3 The accretion disc 

We model the accretion disc as a symmetric, flattened disc 
with user-defined inner and outer radii. In the case of IP Peg 
we forced the inner radius to match the radius of the white 
dwarf. When fitting the individual lightcurves the outer disc 
radius Rdisc was set to match the radius of the bright-spot 
(Rspot, next section) although in the case of the phase- folded 
lightcurves we found a small improvement in the fit when 
we allowed Rdisc to be a free parameter. These fits show 
-Rspot to be somewhat smaller than Rdisc, particularly in 
the tt'-band. This may in part be due to the disc being non- 
circular, but note also Rdisc is rather poorly constrained and 
the uncertainties in these values are likely much greater than 
the formal errors listed in Table [ATI 

In between the user-defined inner and outer radii, the 
height of the disc was defined to follow a power-law in radius, 
i.e. 

h(R) = h R a , (A2) 

where ho and a are parameters of the model. The height pa- 
rameter ho and the radius R, along with all other lengths, 
were scaled in terms of the binary separation, a. We fixed 
ho = 0.02 and a = 1.5. Since the outer disc radius ~ 0.3a, 
this is equivalent to a very thin disc. The feeble emission 
from the disc meant that ho and a have little overall im- 
pact, and so they were not optimised during the fits. The 
surface brightness of the disc was also a power-law in ra- 
dius, with the overall level defined by the temperature at 
the outermost radius of the disc, translated into a surface 
brightness assuming a blackbody spectrum given the cen- 
tral wavelength of the filter concerned. It is important to 
realise that this does not mean that we are saying that the 
radiation from the disc is that of a blackbody; it is just one 
way of normalising the overall contribution from the disc. 
Nevertheless, our model can certainly be criticised for the 
assumption of a power-law variation in surface brightness 
and for its symmetry. Neither is necessarily correct, on the 
other hand there is nothing in the data that definitely points 
the way to a better alternative. We did allow the power-law 
exponent of the surface brightness to be optimised, since it 
has no clear value ab initio for quiescent discs, and since it 
has a significant effect, particularly during the egress phases 
of the eclipse. 

A4 The bright-spot 

The bright-spot was the hardest component to model, and 
indeed, as Figure [2] shows, we never achieved an entirely sat- 
isfactory fit. The bright-spot is of course a uniquely strong 
feature of IP Peg, and so it is to be expected that it would 
cause difficulties. We model the bright-spot as a series of el- 
ements which lie along a straight line in the orbital plane (as 
we illustrate in the Figure schematic). The angle (<j>) that 
the line makes to the line of centres between the two stars 
is a key defining feature of the bright-spot since it is this, 
in conjunction with the distance of the bright-spot from the 
white dwarf and the mass ratio and inclination, which de- 
termines the relative lengths of the bright-spot ingress and 
egress features, which are both well-resolved in our data. 
Each element was then assumed to contribute some fraction 
(fc) that was constant in phase together with some fraction 
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of light modulated according to the angle of the element. 
In our model the elements were assumed oriented perpen- 
dicular to the orbital plane (although we can allow them 
to have arbitary orientation), so they give light-curves that 
peak when we see them face-on. The angle ip is denned as the 
angle away from the perpendicular to the line of bright spot 
elements at which the light from the bright spot is beamed. 
In this manner, the variation of the orbital "hump" and the 
relative size of the bright-spot ingress and egress features can 
be fitted. Note that due to the way they are defined there 
is a strong correlation between the angles <f) and ip. This de- 
generacy can be broken if the observed phase range covers 
the entire orbital "hump" , but this is not the case for our IP 
Peg data. These angles are therefore poorly constrained by 
our model and we see in Table IAT1 that the fitted values for 
(f> and ip vary significantly between the different bands. Note 
however that (fi+tp is consistent, and so the direction of peak 
emission is the same for all photometric bands. Finally, we 
should point out that we find no correlation between these 
parameters and q or i, and so the fact that these bright spot 
angles are poorly constrained does not introduce any addi- 
tional uncertainty in the physical parameters listed in Table 
□1 

The surface brightness of the elements in the bright- 
spot were specified by 



Soc(f)"exp [-(f) 



(A3) 



where x is the distance along the line defining the bright- 
spot, /? and 7 are two more power-law exponents which allow 
some flexibility in how the brightness varies, and I is a scale- 
length. The final element of the bright-spot is its position 
within the binary. The maximum of the surface brightness 
occurs at 



1/7 



(A4) 



We define this location to lie upon the ballistic gas-stream, 
at a specified distance from the white dwarf R spo t, this 
"bright-spot radius" being the final parameter needed to de- 
fine the bright-spot. The path of t he gas stream is calcu lated 
by integration of the equations of iLubow fc Shul l|l975f ). 

Compared to the modelling of the bright-spot in 
iHorne et al.l (| 19941 ') and the related papers mentioned before, 
the parameter 7 allows increased flexibility in the specifica- 
tion of the bright-spot, since in the earlier work, 7 = 1 by 
default. This extra flexibility comes at the cost of increased 
degeneracy; the MCMC method was essential to cope with 
this. Despite this, as mentioned earlier, the strong bright- 
spot of IP Peg is still not modelled perfectly by our code, 
but we felt that adding yet more complexity to the model 
would not be justified by any increased confidence in the 
fundamental parameters such as q and i. 



Table Al. Results from our MCMC fits to the phase-folded 
lightcurves, as detailed in Section 14.11 q: the mass ratio, i: the 
binary inclination. Hi: the white dwarf radius. Rdisc'- the accre- 
tion disc radius. R sp ot- the distance between the maximum of 
the surface brightness of the bright spot (Equation IA4D and the 
white dwarf, a: the binary separation. /3 and 7: the power-law 
exponents for the bright spot. 5: the exponent of surface bright- 
ness over the accretion disc. I: the bright spot scale-length. f c : 
the fraction of the bright spot taken to be equally visible at all 
phases. </>: the angle made by the line of elements of the bright 
spot, measured in the direction of binary motion from the line 
from accretor to donor, tf): the angle away from the perpendicular 
to the line of elements at which the light from the bright spot is 
beamed, measured in the same way as </>. 





tt'-band 


(/-band 


r'-band 


q 


0.457 ±0.005 


0.481 ± 0.003 


0.497 ±0.005 


i 


84.08 ±0.20 


83.64 ±0.10 


83.31 ±0.16 


Ri/a 


0.0058 ±0.0004 


0.0045 ±0.0002 


0.0051 ±0.0003 


Rdiscf & 


0.333 ±0.004 


0.285 ±0.006 


0.281 ±0.004 


Rspot/ & 


0.263 ±0.001 


0.261 ±0.001 


0.257 ±0.0001 


P 


0.51 ±0.06 


0.19 ±0.02 


0.25 ±0.03 


7 


0.727 ±0.022 


0.532 ± 0.006 


0.505 ±0.009 


6 


0.041 ±0.069 


-0.262 ±0.065 


-0.100 ±0.045 


I 


0.0102 ±0.0012 


0.0097 ±0.0001 


0.0096 ±0.0001 


fc 


0.711 ±0.008 


0.723 ±0.004 


0.723 ±0.006 


4> 


105.9 ±1.6 


126.6 ±0.5 


125.1 ±0.6 




59.1 ±1.7 


32.7±0.8 


39.6±0.8 
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